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ABSTRACT 

We present results from high-resolution cosmological hydrodynamical simulations of a Milky-Way- sized halo, 
aimed at studying the effect of feedback on the nature of gas accretion. Simulations include a model of inter- 
stellar medium and star formation, in which SN explosions provide effective thermal feedback. We distinguish 
between gas accretion onto the halo, which occurs when gas particles cross the halo virial radius, and gas 
accretion onto the central galaxy, which takes place when gas particles cross the inner one-tenth of the virial 
radius. Gas particles can be accreted through three different channels, depending on the maximum temperature 
value, Tmax, reached during the particles' past evolution: a cold channel for Tmax < 2.5 x 10^ K, a hot one 
for T > lO^K, and a warm one for intermediate values of Tmax- We find that the warm channel is at least as 
important as the cold one for gas accretion onto the central galaxy. This result is at variance with previous 
findings that the cold mode dominates gas accretion at high redshift. We ascribe this difference to the different 
supernova feedback scheme implemented in our simulations. While results presented so far in the literature are 
based on uneffective SN thermal feedback schemes and/or the presence of a kinetic feedback, our simulations 
include only effective thermal feedback. We argue that observational detections of a warm accretion mode in 
the high-redshift circum-galactic medium would provide useful constraints on the nature of the feedback that 
regulates star formation in galaxies. 

Subject headings: galaxies: formation; galaxies: evolution; methods: numerical 



1. INTRODUCTION 

In the current standard picture for structure formation, 
galaxies are believed to form from condensation of gas within 
the potential wells of dark matter halos (White & Rees 1978). 
Therefore, understanding how gas is accreted onto dark mat- 
ter haloes and galaxies is crucial to understand the process of 
galaxy formation. 

Substantial numerical work recently focused on the 'cold 
accretion' mode (Birnboim & Dekel 2003; Keres et al. 2005) 
as the main driver for gas accretion at high redshift and for 
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haloes less massive than a few times 10^^ M©. The distinction 
between two modes of accretion, cold and hot, was clearly un- 
derstood when the first hierarchical galaxy formation models 
were presented (Rees & Ostriker 1977; Binney 1977; see also 
Benson & Bower 201 1). If the halo virial temperature is larger 
than the temperature of the accreting gas, this will accrete 
supersonically, which causes the formation of an accretion 
shock. If the gas cooling time is longer than the dynamical 
time of the halo, the shock will occur at a radius that is com- 
parable or slightly larger than the virial radius (Bertschinger 
1985; Evrard 1990). If the cooling time is much shorter than 
the halo dynamical time, the shock will occur much closer to 
the forming galaxy. The gas is still heated to very high tem- 
perature, but it cools so rapidly that it cannot maintain the 
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pressure needed to support a quasi-static hot atmosphere . 

This picture has been vaUdated by ID hydrodynamical sim- 
ulations (Birnboim & Dekel 2003; see also Forcada-Miro & 
White 1997), and by recent 3D hydrodynamical simulations 
(Keres et al. 2005; Ocvirk et al. 2008). These simulations 
have shown that a static atmosphere of hot, virialized gas de- 
velops in haloes above 2.3 x 10^ ^M©, and that this 'transition 
mass' is nearly constant as a function of redshift. Simulations 
have also pointed out that, even when a shock is present, cold 
gas accretion can occur along filaments that penetrate deep 
inside the hot halo and might funnel gas to the central galaxy, 
elevating its star formation rates. 

Brooks et al. (2009) investigated gas accretion in a set 
of high resolution simulations spanning about two orders of 
magnitude in halo mass. They find that, for galaxies in the 
range 0.01 to L^, the early growth of the stellar disk is dom- 
inated by gas that is not shocked to the virial temperature of 
the dark matter haloes, and that cold accretion dominates at all 
times for galaxies less massive than our Milky Way. A high 
fraction of gas accretion onto halos is found to be smooth, 
even for the most massive galaxy considered. In a more recent 
study, van de Voort et al. (2011) used the Overwhelmingly 
Large Simulations (OWLS; Schaye et al. 2010) to investi- 
gate accretion rates onto haloes and galaxies, and their depen- 
dence on halo mass and redshift. They confirm that hot mode 
accretion dominates the growth of massive haloes, but show 
that accretion onto galaxies is very sensitive to feedback pro- 
cesses and metal line cooling. In addition, they confirm that 
gas accretion on the halo is smooth. Finally, Oppenheimer 
et al (2010) studied how gas accretion is affected by kinetic 
feedback. They found that a significant amount of the gas 
that later form stars has been previously processed through 
the winds. This wind channel is relevant below z ~ 2. 

In this paper, we re-address the issue of gas accretion 
onto the halo and onto the galaxy, taking advantage of high- 
resolution hydrodynamical simulations of galaxy- size haloes. 
We use our new star-formation and feedback algorithm, 
MUPPI (Murante et al 2010), that includes an effective ther- 
mal feedback but no kinetic feedback. Our aim is to verify 
whether, in this scheme, gas accretion onto the central part 
of the halo remains mainly cold and smooth, and whether 
these results are robust against increasing numerical resolu- 
tion, well above that adopted in cosmological simulations. 

2. NUMERICAL METHODS 
2.1. Simulations 

We study the accretion of gas onto the halo and on its 
central part, the galaxy, using hydrodynamical simulations 
of galaxy- size haloes. Initial conditions were generated by 
Stoehr et al. (2002) adopting a ACDM cosmology, with Q.^ = 
0.3, Qa = 0.7, Ho = 70 km s'^ Mpc'^ and = 0.9. A 'Milky 
Way' halo was selected from an intermediate-resolution simu- 
lation as a relatively isolated halo which suffered its last major 
merger at z > 2, and with approximately the correct peak ro- 
tation velocity at z = 0. Our simulations include as many gas 
particles as high-resolution DM particles, and adopt a baryon 
fraction of ft, = 0.19. In Table [T] we list the main numerical 
parameters of the simulations. The highest resolution simula- 
tion (R2) was run down to z = 2.5. 

Our simulations were carried out using the Tree+SPH code 
GADGET-3, an improved version of the GADGET-2 code 

^ In recent work, gas in cold flows is also coUimated, i.e. it is accreted 
along cosmic filaments 



Table 1 

Simulations. Column 1: Simulation name. Column 2: 
Mass of the dark matter particle, in M©. Column 3: 
Mass of the gas particle, in h~^ M©. Column 4: 
Plummer-Equivalent softening length at z = 0, IN H"^ 
KPC. The softening is physical down to redshift z = 2, 
and comoving afterwards. column 5: number of dm 

and gas particles in the high-resolution region. 
Column 6: Virial mass of the halo, in h~^ M©. Here we 
define the virial mass as that contained in a radius 

enclosing a mean overdensity of 200 times the 
critical density. the simulation r2 was run down to 

Z = 2.5 ONLY. 



Name 


Mdm 




A^DM & Ncas Myir 


Rl 


1.6 X 10^ 


3.0 X 10^ 0.325 


5,953,033 2.70-1012 


R2 


1.7 X 10^ 


3.2 X 10^ 0.155 


55,564,205 



(Springel 2005), by using the star f ormation and feedback 
scheme MUPPI (iMurante et al 11201 Ol MIO). In this scheme, 
each star-forming gas particle is described as a multi-phase 
particle, where cold and hot gas co-exist in pressure equi- 
librium. Mass and energy fluxes between different phases 
are described by a system of ordinary differential equations. 
The molecular fraction of c old gas is estimated using a phe- 
nomenological prescription (|Blitz & Ros olovski"2006). Stars 
explode and inject a fraction / = 0.3 of the available SNe en- 
ergy as thermal energy in the hot and tenuous gas phase of 
neighbours. These do not immediately radiate it away be- 
cause of the hot gas low density. As a result, our scheme 
provides effective thermal feedback: gas particles inside and 
near star-forming regions are heated, and escape from their 
parent structure along the path of minimal resistance (i.e. the 
direction of lowest density), with velocities up to ~ 100 km 
s~\ giving rise to galactic fountains. We stress that we do 
not include kinetic feedback. A particle remains in the multi- 
phase stage for twice the gravitational dynamical time it had 
when 95% of its cold phase was accumulated. This simulates 
the formation and disruption of star-forming GMC. We re- 
fer to MIO for a detailed description of our scheme and for a 
discussion of the choice of the relevant parameters, which we 
kept fixed here. Our model has been included in a recent code- 
comparison project on galaxy formation simulations (Scanna- 
pieco et al, 2011). 

Our simulations also include a phenomenological descrip- 
tion of the UV background (Haardt & Madau 1996), but no 
chemical evolution of the gas: gas cooling is computed as- 
suming primordial chemical composition. 

2.2. Analysis method 

We analyze the thermal history of gas particle accreting 
onto the halo, and onto the galaxy by using the following def- 
initions. The halo is defined as a sphere centered on the most 
bound particle, and with virial radius enclosing an overdensity 
of 200 times the critical density at the redshift of interest. The 
galaxy includes all particles within one tenth of the virial ra- 
dius. We will see that this distinction is very important, since 
the gas accreting onto these components can have different 
properties. To perform our analysis, we used 131 snapshots, 
roughly equispaced in redshift from z = 8.3. 

At each snapshot, we selected the most massive progenitor 
of our z = halo, by using a standard friends-of-friends algo- 
rithm with linking length b = 0.l6. We then assumed that a gas 
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particle is accreted onto the halo (galaxy) at a snapshot n if it 
is part of the halo (galaxy) at the corresponding snapshot but 
not at the previous one n-\. Each gas particle is then traced 
back in time, and the highest temperature reached during its 
past evolution is recorded. 

We distinguish gas accreting in three channels: ( i) cold gas, 
made of particles that never reached temperatures larger than 
^min = 2.5 X 10^ K; (ii) warm gas, composed of particles that 
were hotter than but colder than Tmax = 10^ K; (Hi) hot 
gas, made of particles that reached temperatures larger than 
Tmax- In previous work, the warm channel has been included 
in the hot one. In this study, we have considered this separate 
channel to quantify the accretion that occurs at 'intermediate' 
temperatures, that are indicative of feedback from supernovae 
(SNe) rather than shock heating. 

Finally, we identified all self-bound structures at each snap- 
shot using the algorithm SUBFIND (Springel et al. 2001), and 
built merger trees for all of them using the same procedure 
described in Springel et al. (2005). 

Using these merger trees, we check if gas particles accreted 
on the halo or on the galaxy belong to a distinct subhalo at the 
previous time, i.e. for each accretion channel, we determine 
the fraction of gas accreted in "clumps". Gas particle that 
were not assigned to a subhalo in the past are considered as 
"smoothly" accreted. 

3. RESULTS 

Figure [T] shows the fractional mass accretion rates in the 
three channels defined above. Right panels are for accre- 
tion on the halo, while left panels show the corresponding 
accretion rates on the galaxy. We show the accretion rates 
for the resolution levels Rl and R2, as well as for a Rl run 
that adopted the Springel & Hernquist (2003) effective model 
(see below). The Figure shows that gas is accreted on the halo 
mainly via the cold channel for all runs considered. At low 
redshift, z < 2, the importance of warm and hot channels in- 
creases, in agreement with previous studies (as in, e.g., Keres 
et al. 2006, ,Brooks et al 2009, van de Voort et al 201 1). 

Results are quite different for the accretion onto the galaxy. 
In this case, the cold and the warm channels equally con- 
tribute to gas accretion at high redshift. At later times, the 
hot channel provides the dominant contribution. 

Our findings are robust against resolution: our R2 simula- 
tion has a mass resolution which is sufficient to resolve halos 
of ^ 10^ M0 with ^ 100 DM particles. Therefore, results 
from this simulations are reliable also at high redshift and, as 
shown in Figure [B they are in excellent agreement with those 
obtained for the Rl run. 

Figure[2] shows density maps of gas which accreted onto the 
galaxy between redshift z = 2.5 and z = 2.0. We show the gas 
at previous redshifts: z = 3.8 (left column), z = 2.9 (middle 
column), and z = 2.5 (right column), for the Rl simulation. 
Maps are centered on the most massive progenitor at the cor- 
responding redshift. The upper, middle and lower rows are 
for the cold, warm and hot channels, respectively. The Fig- 
ure shows that the cold gas is characterized by a very clumpy 
distribution, while the hot gas is always diffused around the 
most massive progenitor. Interestingly, warm gas is more dif- 
fuse than cold gas, but it surrounds the cold clumps. This sug- 
gests that the infalling cold clumps are strongly star forming, 
and that they heat and eject part of their gas while infalling 
towards the central galaxy. 

This is confirmed by the analysis of the multi-phase nature 
of the accreting gas: SN thermal feedback is effective at heat- 



ing multi-phase particles. So, by computing the fraction of 
gas particles that were multi-phase when they reached their 
maximum temperature, one can disentangle the influence of 
SN heating from gravitational shock heating. For the Rl sim- 
ulation, averaging over the entire redshift range (from z = 8.3 
to z = 0), we find that when accreting gas reaches its maximum 
temperature, ~ 97% of it is multi-phase and star forming in 
the cold channel. The corresponding fractions for the warm 
and hot channels are ~ 66% and ~ 8%, respectively. This 
suggests that gas in the warm channel is heated by SN feed- 
back, while gas in the hot channel is heated by gravitational 
shocks. 

To better illustrate the role of effective thermal feedback 
we run simulation Rl using the widely used effective model 
of Springel & Hernquist (SH, 2003) for star formation, with 
kinetic feedback and velocity = 340 km s~^ ^. The result- 
ing accretion rates are shown in the lower panels of Figure 
[TJ the contribution of the warm channel to accretion on the 
galaxy becomes negligible, while gas accretion in the cold 
channel increases. This confirms that simulations with unef- 
fective thermal feedback are unable to heat cold flows into the 
warm channel, even when they produce significant outflows. 

Finally, for each channel, we analyse the dumpiness of gas 
accretion, both on the halo and on the galaxy by determin- 
ing if its particles belonged to an accreting satellite of were 
accreted smoothly. 

Results are shown in Figure [51 For accretion on the halo, 
we find that total gas is mainly smooth (as in, e.g. Brooks et 
al. 2003, van de Voort et al. 201 1). At z = 0, 69 per cent of the 
gas has been accreted smoothly in the Rl simulation. Smooth 
accretion dominates the hot channel, while in the cold one, 
smooth and clumpy accretion contribute equally to the total 
accretion rate: at z = smooth accretion accounts for only 49 
per cent of the overall cold accretion. 

Also the total accretion rate on the galaxy is mainly smooth, 
^ 67 per cent for Rl, and hot accretion shows the same trend 
found for the accretion on the halos. We find a significant 
difference, however, in the cold channel. In fact, this channel 
is dominated by clumpy accretion, with 60 per cent of the 
gas accreted in this way. The gas accreted through the warm 
channel is also prevalently clumpy. 

We stress that the clumps that contribute to the cold and 
warm channels have a significant (> 85 per cent) DM frac- 
tion, ruling out possible numerical problems of the SPH tech- 
nique that can, in certain conditions, produce "cold-gas only" 
spurious blobs. 

4. DISCUSSION AND CONCLUSIONS 

We analyzed the gas accretion rate in cosmological numer- 
ical simulations of a galaxy- size halo. We take advantage of 
our new MUPPI star formation and feedback model that pro- 
vides effective thermal feedback and use different resolution 
levels to test the robustness of our results. 

We distinguished between gas accretion onto the halo (de- 
fined as all the matter included in the redshift-dependant virial 
radius), and accretion onto the central galaxy (defined as all 
the material included within one tenth of the virial radius). 
We traced back all gas particles in our simulations, and de- 
fined a gas particle to be accreted at redshift Zn if it was not 
part of the main progenitor (or of its central galaxy) of our 
z = halo at z«-i > z«, but became part of it at z«. We defined 

^ We decouple wind particles from ambient gas for 20 kpc or until density 
drops to one half of SF threshold 
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Fig. 1 . — Right panels: gas fractional accretion rates onto the halo, as a function of the accretion redshift. Left panels: corresponding fractional gas accretion 
rates onto the galaxy. In red, we show the hot gas accretion, in green the warm accretion, and in blue the cold accretion. In pink, we show the total gas accretion 
rate. For clarity, we have resampled our redshift outputs in 20 redshift bins. 



three accretion channels: a cold one, including gas particles 
whose maximum past temperature was lower than 2.5 x 10^ 
K, a hot channel, that includes gas particles with maximum 
past temperature larger than lO^K, and an intermediate warm 
channel. 

Our results can be summarized as follows: 

• The warm accretion channel becomes as important as 
the cold one at high redshift, when considering the 
accretion on the central galaxy. This trend is robust 
against numerical resolution. 

• Gas particles in the warm channel are heated by SN 
feedback. This is confirmed by the analysis of the 
multi-phase nature of the gas, and by the fact that a 
simulation of the same halo carried out with a model 
that does not include effective thermal feedback shows 
negligible accretion rates through this channel. 

Regarding the accretion on the halo, our findings confirm 
previous studies: the 'cold' accretion mode on the halo is im- 
portant at high redshift, while the 'hot' accretion mode be- 
comes dominant at lower redshift. However, accretion rates 
onto the halo are different from accretion rates onto the central 
galaxy. We find that a significant amount of gas is accreted on 
the galaxy through a warm channel. The gas is in cold clumps 
when entering the halo, and becomes warm because of the SN 
feedback originating from star formation within the clumps. 

Van de Voort et al. (2011) too analyzed the dumpiness of 
the accreted gas, but they included all clumps with mass ratio 
smaller than 1:10 in their smooth component. These clumps 



are, however, very dense and star-forming, and effectively 
heat the cold gas to the warm channel. 

Oppenheimer et al. (2010) found that a large fraction of 
low-z star formation experienced a wind phase in the past. 
However, they studied the effect of kinetic feedback, and 
showed that this channel becomes important below z = 2. 
Their results are thus complementary to those shown in this 
letter. 

The picture emerging from our analysis can be sketched as 
follow: at high redshift, accreting gas is mainly cold, partly 
clumpy and partly smooth. Not all of the smooth cold gas 
reaches the central part of the halos because, behaving as 
non-collisional material, cold, un-shocked gas is not able to 
loose orbital energy and angular momentum. Thus, only gas 
streams with very low impact parameters can funnel gas to the 
central galaxy, contributing to its star formation. Cold clumps 
suffer dynamical friction and spiral towards the center. Since 
they are very dense, their star formation rate is high, and part 
of the gas they contain is heated by SN feedback, and is ac- 
creted onto the galaxy as warm, slightly less clumpy, mate- 
rial. We do not ascribe these findings to specific features of 
our sub-resolution model, and argue that any effective thermal 
feedback scheme would give similar results. In this picture, 
the time needed by the cold gas to reach the central galaxy is 
neither the free fall time nor the cooling time. Smooth cold 
accreting gas can only reach the galaxy if the impact parame- 
ter is small. In contrast, cold clumps are subject to dynamical 
friction so they can reach the central galaxy even with non 
negligible impact parameters, over a time- scale that is deter- 
mined by dynamical friction. 

However, the exact quantification of warm flows may well 
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Fig. 2. — Projected density maps of gas accreting on the galaxy between z = 2.5 and z = 2, for the Rl simulation. Different rows correspond to the three 
channels considered: cold (upper row), warm (middle row), and hot (lower row). Different columns correspond to different redshift: z = 3.8 (left column), z = 2.9 
(middle column), and z = 2.5 (right column). The scale of the maps are 960, 480, and 240 h~^Mpc comoving for the three columns. The color scale is logarithmic 
from 1.0 to 10^ Mq h pc^ comoving (from pale pink to red). 



depend on model details. In particular, the entrainment of 
cold gas into the multi-phase particle may lead to some over- 
estimation of the amount of warm flows. We will address this 
issue in future work. 

It is widely accepted that stellar feedback plays an impor- 
tant role in determining the physical properties of galaxies. It 
remains, however, unclear if the nature of feedback is predom- 
inantly kinetic or thermal. Our results demonstrate that stellar 
feedback strongly affects the termodynamical state of the gas 
accreting onto galaxies. In particular, we argue that a mode of 
feedback that is primarily thermal would produce significant 
amount of warm gas accreting onto galaxies at high redshift. 
This has noticeable consequences on the observability of gas 
flows infalling onto high redshift galaxies, and may explain, 
at least in part, why cold flows have been difficult to detect, to 
date. 
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Fig. 3. — Smooth vs clumpy accretion. Solid lines refer to accretion on the galaxy, while dashed lines are for accretion on the halo. Blue and red lines 
correspond to clumpy and smooth accretion, respectively. The upper panel shows the total accretion, while the three channels considered are shown in the lower 
panels: from the 2"^ upper to lower panel, accretions for hot, warm and cold channels. Similar results are found for lower resolutions (see text). 
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